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STRACT 


There is much disagreement as to (a) what causes clear air 
turbulence (turbuience which is not in or near convective clouds and 
is above 15,00) feet in altitude) and (b) which seocercloetea param- 
et2%s tan ke used to detect and forech its occurrence. The approach 
to this problem has beer to relate not one parameter to clear air 
turbulence but various parameters, By summing these parameters areas 
can be defined where there is a high probability of encountering clear 
air turbulence. Each parameter has been based on a statistical study 
which found a relationship with clear air turbulence. The parameters 
used were horizontal and vertical shear, curvature, kinetic energy 
and their derivatives. | 

The numerical forecasting program proposed ere can be extended 
to the stratosphere when more reliable height and temperature fields 
are available. This program will have much more significance when 


intermediate forecast height fields, temperature fields and a grid of 


much smaller mesh length are available. 
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sue UPTRODUCT LOB 


The phenosenon of clear sir turbulence (henceforth denoted as CAT) 
appesrs to consist of random three dimensional eddies that occur in 
certain confined regions of the *tmos3phere. This phenomenon was firet 
Tecegnized in the early 1940's as “air pockete'', and has gained in im- 
portance ax depth of study with the development of fast-flying swept 
wing and delta winged alrcraft. It is of prime importance te the avia- 
tion industry vecause it affects the safety and cemfort of passengers 
and crev, as weil es operational costs. There have been several cases 
where loss of coatrol of aircraft, structural damage, passenger and 
crew iniurfes (even dzaths) have resultei from CAT. 

Airersft manufacturers and the airlines are facerested in CAT in 
order to determine the strength of airframe design so that it can be 
built to withstand all degrees of turbulence. Planning of supersonic 
transport aircraft is considering the effects of this unexpected turb- 
ulence occurring at auy flignt level {4n the atmosphere. Therefore, it 
becomes an even greater operational prevlem in the era of future design 
of SST aircraft as to cost, safety and comfort. CAT is usually less 
intense than turbulence encountered in thunderstorms 1 . CAT may be 
considered at times more dangerous than “thunderstorm turbulence" as 
it usually occurs with no visible warring. 

The various military agencies ere deeply concerned with the under- 
standing of CAT and ites prediction so they may be better prepared to 
take it into account operationally in all areas of the world. Scientieta 
are highly interested in CAv because it is a phenomenon of our environ- 
ment that {s not clearly usderstood. We are faced with the problem of 
heving to rely on mesoscale data while the phenomenon is of a microscale. 


Because of the great differences in scales involved in the forecasting 


difficult te identify. It is the intent of this paper to develop a 
numarical forecasting system of CAT whereby an evpirical approach has 


t 
problem, isclated patches of CAT (in eithes space or times) ere vary 
heen used to escomplish this end. 
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2. DATA COLLECTION 


The nature of clear air turbulence, its physics anc its meteorology, 
sre still not completely maerriood. More mesoscale ‘Zim to 10Cm 
hborigomtal distance) amd ricroscele (less than 2is) studies mit be 
mace comcerning (he cauces am! generation of CAT. Any improvement in 
forecasting skill thst say be realized Sy this study wiil prebably result 
from increased mecoscaie input data into the forecest problem. The only 
mesozcale measuremects available opersticnaily at this time which concern 
the problem are rawinsonde cata that give a neariy concinuouz zeasure- 
went in the vertical of wind, presscre, temperature and humidity ise 
the surface to above Wie ir altitude, and satellite clow! pictur2s 
thet deliseate clus strmetures down to a scale of about Jim. 

The uzual detinition of CAT ic "atsospheric turbuience whica is 
not im cr near tonvective clouds, including thuaders*omms, and is net 
below 15,000 feet in altitude’, ‘4) Thus, mechanical turduienc: induced 
by rough terrain is exciudec. It is realized that thie is an extremely 
srbitrary definition resulting from a desire to simplify pilot reporting 
procedures. 

Turbulence intensities are, at present, designated as light, mod- 


erete, severe end extreme. In recent yez.s severcl cefiaitions of the 
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four cetegories of turbulence have been proposed. The current defin- 
itions in officic! use were develcped by the National Advisory Ccmmittee 
for Aeronautics (RACA) Subcommittee on Meteorological Froblems (1957). 


See Table 1. 


There are few quantitative measurements of atmospheric turbulence 
at any scale. The buik of the aircraft turbulence data available is 


based upon the four categories stated in Table 1, (4) , and, conse- 


quently is highly subjective and qualitative in nature. Included in the 


factors that affect the sudjective decision of the pilot when reporting _ 


kurbulence are the wing loading, the aircraft's speed and attitude, ehe 
pilot's training, experience, and his ey choi oeitel Peetiens Ce. 
Objective criteria for turbulence reporting must vate until more 
Measurements are made of the conditions existing in tne turbulent eddies, 
and watil cockpit instrumentation includes a fet load a turbulence | 
indicator. 
Much time, eftort and finar.cial support have been expemed in thie 


country by the govermment, the aircraft industry, and private ond public 


research facilities to investigate the problem of clear air turvulence . 


and to find a good approach leading to a celution to this prublem. This 
research has gone forward in three broad categories: first. eee 
concerning the correlation of CAT with macrogcale, mesoscale, and/or 
microscale atmosyheric measurements; secend, research into inst rumea- 
tation fcr the detection of CAT sufficiently far in advance of the air- 
craft to allow evasive measures (6) and third, research by military and 
airline organizations concerning tne operaticnal aspects of CAT. 

Various projects such as HICAT, ALLCAT, and TOPCAT Have een unde =- 
taken to study clear air turbulence. The problems encountered were many 
but the results showed that it was indeed possible, witt suitable in- 
strumentation, to find, track and record CAT. The use oi cocetant level 


bsllons and doppler radar to detect CAT nas met with some st:cceas, but 


again tie lacx of detail necessary in CAT studies leaves imch to be 


desired. 
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Measurements of the microstructure, which contains the perturbetions 
of CaT dimensions (100 to 569 meter wavelength), have been taken by 
specially inetrumented aircraft. Until recentiv the aircreft itseif 
was used a8 a sensor to sseasure atmospheric gusta from the aircraft 

i acceleration dats. The Air Force realized the danger of relying eclely 
ou the acceleration response of the U-2 aircraft as @ measure of turbu- 
lence at high altitudes. A full knowledge of che aircraft's response 
to turbulence over a wide range of wavelengths is required for maaning- 
ful imterpretation of such data. This procedure usuaLly gives reliable 
results at short warclengths up to a few hundred metere. At long wave- 
lengths, this instrumentation becomes ies sensitive. Accclerations in 
Longer waves are usually small and can be masked by pilot induced air- 
craft motions. 

4 .| & supersonic or hypersonic craft of some radical shape, flying four 

to ten times the speed of the U-2, will obviously have a somewhat dif-~ 
ferent: respanse to turbulence than the U-2. An aircraft flying at 

shese hig.: speeds would be affected much more by longer turbulence wave- 

lengths aad ieas by the shorter. 

As pointed out earlier, the meteorologist has few direct measvie- 
ments of turbulence intensities and must depend upon the - -curacy of 
the intensities reported by pilots. In our study, we used the Colson 
monthly turbulence summaries which were obtained from the Air Force 
3rd Weather Wing at Offutt AFB, Omaha, Neb. This rep..t suited our 
needs most closely because the CAT sceports were detailed as to location, 
time, altitude and intensity, In the period used in our study, Si:cember 


1964 through March 1965, there were no less than 3670 CAT reports. The 
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3. TRE CLEAR AIR TUKBILERCE STUDY 

Clesr air turbulence (CAT) is a microscale phenomenon (less than 
2 kilometers in horizontal extent) but the conditions which are aympto- 
matic of .t- exieteace are of anyopiic scale. Therefore, synoptic 
seale parameters =ay be uzed to datermine areas where CAT could occur, 
that is forscessting areas where there 1e a high probability of CAT. 

From various reporta and thei: contradicaicns it appears that no 
sircele parameter can detect CAT. Certain parameters can detect the 
possisility of CAT in some syneptic situations but fail in others. 

The approach used in this paper was to take statistical studies 
made by previous investigators and to use the basic parameters which 
they related to CAT (2,3,5,7,8). Lf the magnitude of any one parameter 
becomes large or the sum of various parameters becomes large then there 
ig a high probability of CAT in that area. Therefore, the problem is 
net one of forecasting actual CAT but rather to forecast areas of high 
and low probability of encountering CAT. In this way flights may be 
planned so as to expect least CAT. 

The research program consisted of computing the equations shown 
in Appendix A. These equations were computed using the CDC 1604 digital 
commuter. The program was written utilizing symbciic coded relocatable 
assembly program (SCRAP). It was necessary to use fixed point fractional 
muabere in order that Fleet Numerical Weather Facility (FNWF) subroutines 
could be used. Ail finite differences were computed using standard FNWF 
mesh length of 381 kilometers true at 60 degrees north latitude. There 
ere no time derivatives in the program. The research reported here is 


accomplished using analytical fields, Operational use would employ 
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forecast fieids. Results would naturally net have been as good if fore- 
cast fields were used in thig research. ‘Prog 24 hours" printed st the 
bottom of each field represents the practical forecast interval to be 
used operationelly. 

The research program was written to compute CaI in turee Larers 
500 to 300, 300 to 200, and 200 to 100 miilibars. Computations were 
not made for the third layer since 100 millibar fields were not available. 

The research program was written to compute on the entire €3 x 63 
FNWF grid of the northern hemisphere. A boundary condition of zero was 
used for the outside vows and columns. The print routines are 22 x 22 
extracts of the United States starting at FEWF grid point J008, 1016. 
The latitude and longitude coordinates of the four corners of the printed 
fieids starting with the lower left corner proceeding clockwise are: 
9.1N 109.5W, 44.5N 165.5W, 57.78 3.18, 13.0% 58.8%. 

Ali printed fieids are pure mmbere and have so disersional 
meaning. All scaled outputs have been chifted to the left end of the 
register and the first three mmbers with sign bit are printed oat is 
decimal. Grid points can take on values from -999 to +399 except those 
fields which have been made all positive. The decimal point does act 
appear on the printed fields, Thezefore, the printed grid point values 
are from -999 to +999. 

In the research progras: the capital letters refer to the paramster 
as computed from the data. The small case letters serving es expozents 
represent the number of times and direction t’¢ register has to be 
shifted in order to place the significant portion into the first three 
numbers. Therefore, the two with exponents represent the scaling co- 


effictent. Since the computations were in fixed point fractional e221 
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printests aed to be shifted s> as Bot to exceed ome st any srid soint 
in che field. @zxeeeding ope wocld revzir in 4 meaningizes vakue ast 

thet parciculer grid point. [ea addition the entire rieid bad to be 

kere large escugh so that patteres conit exist acd 20% be st or pear 
sere threwgheut the fielé. Since toese fields are summed they west be 
omal] omomg> so 26 sot to carse che summation field ro exceed ome at 
exp grid pcine . 

The rerearch progres was sem for thircy-siz days dcricg December 
LPS ami lemeary, Febreary, end Merch i865. Tee thirty-six <days were 
chosez. becsase Chey were Se aost active is CLAY reports éuring tee four 
mom potied. fe ether words, there were morte reports br pilots ex- 
comexevrimg CAI ce theee paerticeiax: <eys. is order to show 22 exemic 
ef each Held pristet oat ty the recearch program the tise OOl 23 FEB 65 
wae surbizrarily selected. Theme fieids appear ic Appecrtix ¢. Eaca field 


was profeced by ac imtividua! term weict will be cescribed as follows: 
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The Firet Term ts 2* 4 where i is che Leclacian of absciete 
eortisizy. Weer this twze is segetive sbere is 2 iccet mexinam of 
s@soigcte worticicy meaning (¢ is lasger 2t tear gris poist chee che 
everauge of the parimemicg grit poicts. Tetrefore, the creloaic curvwe- 
Tere sr cypciocic stear or bat are relatively large ec that grid soit. 
Teis chocift correrpont <5 tee coif site of che set espectaliy ic 
throughs. Acooeding <5 Entlict axe Melee (1) tiese is a greeter per- 
comtage <£ Cat oc the colt side of the yer. Also acscorting to Rarrisan 
(Z) Cheve is < stroeg cemtemcy fox moderate 25 severe CAT to be assc- 


ciated sieh irceagt iimes. 
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This term wes computed for the layer by first caiculeting the 
absolute yorticity of the upper end lower level D fieids. The Laplacian 
was then taken of each fieid and a vertical average made of the upper 
and lower levels to obtain the Laplacian of vorticity of the layer. 

This parameter hae sepicted most of the CAT associated with the 
trough over the western United States. The severe CAT near New Orleans 
is in an area of very large negative aumbers. The field 4s contoured 
at intervals of 100 with the origin at zero. 

THE SECOND TEER 

The Second Term is 2° B where B is the absolute value of the 
vertical chenge in the vector thermsl wind. One of the parapetezs 
which Lake's (7) statistical testirg indicated was associated with CAT 
was the vertical gradient of wind shear. As showa by Richardsen (9) 
the thermal wind shear is proportional to the gradient of static 
stability. 

The u and v componenrs of the thermal vind were computed from the 
upper and lowe: ievel temperature fields. The difference between che 
upper ard Lower level values of the u component was found and each, 
difference was equared. This was also done for the v component. The 
square root of the sum of the squared differences piven ee aaentcede 
of the vector difference. According to Endlich and McLean (3) the 
largest values of the thermal, wind shear appear on the warm side of the 
jet. This was found to be true throughout the foux months. The contour 
interval for this field is 25 and the origin is seio, 

THE THIRD TERM 

The Third Term is 2© C where C is one half the geostrophic wind 

velocity squared. C therefore represents the epee tic kinetic energy 


or in other words the kinetic energy per unit maser 


16 


cern 


Clem (2) found that most cases of moderste to severe CAT were 
associated with areas of isotach maxima. 
This ters was computed for the layer by calculating the u and v 
: components of the geostrophic wind at the upper and lower levels. The 
upper and lower level u components were vertically averaged to obtain 
an average u component for the iayer. The average v component for the 
layer was obtained by a similar process. The magnitude of the velacity 
Sinhed sis just the sua of the squared components. 
In the research program this field is contoured at intervals of 
100 with the origin at zero. Contoured at this interval the kinetic 
| energy field clearly depicts the isctach maxima regions. The kinetic 
energy field in Appendix C shows this field depicting the CAT in the 
: western part ef the United States occurring in areas of relatively 
, ’ large wind speeds. The xinetic energy field fails to indicate the 
severe CAT near New Orleans because it occurs in an area of relatively 
-_ Light winds. 
THE FOURTH TERM 
The Fourth Term is 24 D where D is the absolute value of the 
derivative of the kinetic energy with respect to pressure. 
Lake's (7) statistical testing indicated that the gust intensities 
are related to the vertical gradients of horizontal kinetic energy. 
This ete was computed for the layer by first calculating the velocity 
eresred ae tte upper and lower levels. The vertical gradient for the 
layer was then pbeained by computing the difference between the upper 
and lower level values of the velocity squared. The absclute value was 


taken so as to have all values positive. CAT should be assocfated with 
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large values ef this field. Large values of this field were found 
only in areas of large values of kinetic energy. This field was 
therefore redundant and vas eliminated from the CAT forecast program. 


The conteur interval was 100 with the origin at sere. 


THE FIFTH TERM 

The Fifth Term is 2° E where E is the absolute value of the 
Laplacian of kinetic energy. The statistical survey made by Endlich 
and McLean (3) shows the maximum occurrence of CAT along che edges of 
the isotach waxima, The Laplacian of kinetic energy ehows Large 
horizontal changes in kinetic energy, both positive and negative. 
Therefore the absolute value of the term is taken in order to give 
only positive mumbers. The contour interval was 25 with the origin 
at zero. - | 

This term was introduced to depict the areas of iarge horizontal 
Skane in kinetic energy. However, there was no relaticnship with 
reported CAT. This field was therefore eliminated from the CAT fure- 


cast progran. 


THE SIXTH T2RM 

At this point in the research program it was necessavy to aum the 
first five terms due to computer memory space. This term, referred to 
as KATL, was the Sixth Term. The contour interval was 250 with oe 
erigin at zerc. It was found that this term did not eupply signif- 
icaptiy new information since it was dominated by the kinetic eae 


+ 


and the two associated terms, 


THE SEVENTH TERM 
The Seveach Term is 2£ F where F is the Jacobian of temperature 


and omega (component of the wind norwai to the pressure surface). This 
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term wea developec by Dr. Moore of Douglas Aircraft and Dr. Krishnamurti 
(>. Te Latter wee associated st that time with the University of 
Califernis, Los Aageies end cocsultant to Douglas Aircraft. The term 
was cevelopec eas the Jaccbian cf cemperacure and three dimensions! 
éivergence. As shown in their peper this is propertional to the nega- 
tiwe of the Jacobian of temperature and omega. This term was computed 
Yor the rower lewei cf each Layer in the research program. The contour 
feterval for this terw wes i100 with the origin at zero. In this pro- 
Kram cc significesc relsttomship wae found vith Large negative or 
positive sarbers over the fovr month period. This tense was therefore 


elissteszecd from the CAT forecest pregran. 


TE ey Tee 

The Zighth Term is 28 G where G is the absolute value of herizontal 
ciwergesce. It was compoted by taking the derivetive of omega with 
respect to presstre. As previovsly strated the entrance and exit regions 
of isctach mzzima areas have been found to be associated with CAT. 

Yheee areas are aiso associated vith horizontal convergence at the 
entrance end divergence at the exit regions. Therefore the absolute 
valce of the chenge in omega with preasure represents the convergence 
eed divergence in the layer parallel to the pressure surfaces. 

Teis term vas computed by subtracting the lower level omega value 
Ero the copper Level amege value at each grid point. Areas of can- 
vergence mad divergence of the height fleld are quite vividly depicted 
by the diverge.se fleid. Ho significance relaticcskip wae found between 
che aivergesce field? acc the CAT -rcurrences, therefere, it vas elimin- 


atec from the CAT forecast program. 
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THE NINTH TERM 

The Ninth Term was KAT2,the summation of all previous terms. Ro 
significant releticnship was found with CAT occurrences since several 
terms tended to cancel each other out. 

. The program then goes into the second layer from 300 to 200 milli- 
bars. All terms were computed and sealed the same, except the diver- 
gence term. [It could not be computed becauee the 200 millibar omega 
field was rot available. 

The third layer from 200 to 100 millibars could not be run for 
these time perfois since the 100 miilibar fields were not available. 
There were very few CAT reports above 200 millibars, therefore, the 
loss was insignificant. 

In all three layers the lower level height field 3 Stace out 
Fite This gives a general impression of the synoptic situation and 
renders more significance to the patterns developed in the other fields. 
The contour interval for the 500 millibar field is 60 meters with he 
origin at 5580 meters. The contour interval for te 300 millibar field 
is 120 meters with the origin at 9120 meters. The contour ouetral for 
the 290 millibar field is 120 meters with the origin at 11,760 meters. 

The GAT forecast program appears in Appendix D, The first three 
terms of the research program are used with ainor changes. “In the flies 
term "al is changed to minus one and the contour jae val bes Bran 
changed to 150 to give better defined pectoris The second re has 
Recnticeed a The third term is iachantedeereene for the con- 
tour interval which was changed to 50 to incress2 chau pecteen size. | 
The KAT field itself is the sumaation of these three rs <a gives 


quite reasonable patterns and pattern sizes. 
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The pattern area depicting ae high probability of CAT would neces- 


H 
sarily be larger during a more active CAT period. The most active | 
part of the year was the four month period December 1964, January, | 
; February, and March 1965. During these months the most active three 
day period was the tenth through the twelfth of March 1965. Therefore | 
the KAT fields have quite large pattern sizes in Appendix E since they 
repreeent the post active three days of the year. The XAT field 
patterns were smailer for less active periods. The total area covered 
by these patterns is much lesa in the KAT field than in the other thzee 
parameter fields. 
This is exactly what was attempted in order to obtain optimum 
size of the forecasted danger areas. If the KAT field patterns ere 
too.large, flights will be rerouted unnecessarily. On the other hand, 


, 2 if the KAT field patterns are too smail, there is a reai danger of CAT 


occurring outside these areas. Therefore, the restraint of the KAT 


field patterns is necessary in order to have an operaticnally useable 


preduct. 
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4. DISCUSSIGH ANP RECHORIDATION FOR FUTURE STIDIES 

The period used in this paper was chosen decense of the largest 
nusber of reported CAT occurrences. Sf the four month peried (December 
1964, Jamusry, February snd March 1965) there were scattered pericds 
where ea large number of eccurrences were reported. This four-month 
period was used to dstermim which parameters were beet suited for 
forecasting purposes. The three-day period of 10-12 March 1965 wes 
selected for a correlation atudy in order to find out how successful 
our forecast sethod is. Tables 2, 3, and 4 show the various fields 
used and the resultant KAT fleld for the 10th, 11th and 12th 5¢ March 
1965. 

The use of the term "percent correlation" as used in this study 
does not mean to imply a statistical correiation. The ideal forecast 
verification makes use of those cases where CAT is forecast, but does 
not occur, and where CAT is not forecast and does not eccur. Ia our 
study it was impossible to take those cases quantitat‘.vely into accovnt. 
Therefore, it is to be understood that "correlation" as used in this: 
study was a general comparison of those reported CAT occurrences that 
fell within the delineated area of high probability of CAT against those 
that did not. For example, if there were ten reported CAT occurrences 
for a given layer and time reriod and six of these reports fell inside 
or on the line delineating the CAT urea and four reports izil outside 
the area, then for that field, layer and time period we would list it 
as six occurrences correlated or a sixty percent correlation. 

Listed are the names of the fields, the number of occurrences 
of CAT for each field aud the percent correlation by field. Also shown 
are the number of occurrences and the number of occurrences that cor- 


related by field aml CAT intensity. 
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After eombiring the Laplacian of Vorticity, Vertical Cradiem cf 
Therasi Wini and Kinetic Energy we arrive et che KAT Field which is oer 
j end product for the area of high probability of CAT occurrences. Zves 
though « higher correlation may be seen in some fields other than the 
KAT field, cwe gbould realize that these fieids encompassac a larger 
than sverags area. In such ceses cue must expect a high correlation. 
The high correlation 12 then not due to the finesse of the forecast 

method, but rather due to the fact that for most af the DSA these wac 
& forecest of high provsbility for CAT. Theareticaily, it would he a 
good idea to divide the percent correlation (such as we computed} Ly 
the size of the area for which CAT wae forecasted. Ye did not fallow 
this idea quantitatively, but caly qualitatively. Therefore, ome finds 
thet the percent correlation for max ultimate forecast (labeied EAT} 

’ a8 sometimes lower thac the percevt correlation for owe cf tre three 

separate fcrucastimg fielas. 


Table 5 is s summary for the three day period. It shows the total 


f numer of CAT occurrences by turbulence category, percent correlakien 
‘ ty field and turbulesce category, end the three day percent carreiaticn 
vy field. 


| _ Phe xresulte were most encouraging ani we fecl thar our emi product 
was a substantial step in at least the right direction toward foreersting 
clear air turbulence. Our knowledge of the meso-~ and micro- structure 
of flc catterns in tie free atwosphere, especially above the tropopause, 
is rill rather poor. A strong need still exists for a well organized 
and well equipped azesurement program, especiall~ at fiight levels of 


[ the future supersonic transport aircraft. Measuiement progress using 
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{3 ABSTRACT 
There is ach disagreezent as to (aj what causeg clear air turbulence 
{surbulence which is not in or near convective clouds and is above 15,000 feet 
in sltitude) and (b) which meteorological parameters can be used to detect and 
i | forecast its occurrence. The approach to this problem has been to relate not 
one perameter to clear air turbulence but various parameters. By summing these 
parameters areas can be defined where there ie a high probability of encountering 
‘ clear air turbulence. Zach parameter has been based on a statistical study which 
: fourd a relationshiy with clear air turbulence. The parameters used were hori- 
* f zontal and vertical shear, curvature, kinetic energy and their derivatives. The 
anumerical forecasting program proposed here can be extended to the stratosphere 
when more reiiable height and temperature fields are available. This program will 
ep have wuch more significance when intermediate forecast height fields, temperature 
# fields and a grid of much smaller mesh length are available. 
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